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Abstract 

Properties of the fractional Schrodinger equation have been stud- 
ied. We have proven the hermiticity of fractional Hamilton operator 
and established the parity conservation law for the fractional quan- 
tum mechanics. As physical applications of the fractional Schrodinger 
equation we have found the energy spectrum for a hydrogen-like atom 
- fractional "Bohr atom" and the energy spectrum of fractional oscil- 
lator in the semiclassical approximation. A new equation for the frac- 
tional probability current density has been developed and discussed. 

We also discuss the relationships between the fractional and the 
standard Schrodinger equations. 

PACS number (s): 03.65.-w, 05.30.-d, 05.40. Fb, 03.65. Db 

Keywords: path integral, Levy motion, fractional differential equa- 
tion. 



1 Introduction 

The Feynman path integral approach to quantum mechanics [|T| , is in fact 
integration over Brownian-like quantum mechanical paths. The Brownian 
motion is a special case of the Levy a-stabie random process. In the mid of 
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1930's P. Levy and A.Y. Khintchine posed the question: When the sum of N 
independent identically distributed quantities X = X\ + X 2 ... + Xn has the 
same probability distribution p N (X) (up to scale factor) as the individual 
steps Pi(Xi), % — 1, ...AT? The traditional answer is that each Pi{Xj) should 
be a Gaussian, because of a central limit theorem. In other words, a sum 
of iV Gaussians is again a Gaussian. P. Levy and A.Y. Khintchine proved 
that there exists the possibility to generalize the central limit theorem ||, 
[|]. They discovered class of non-Gaussian Levy a-stable (stable under sum- 
mation) probability distributions. Each a-stable probability distribution has 
a stability index a often called by the Levy index, < a < 2. When a = 2 
the Levy a-stable distribution is transformed to the well-known Gaussian 
probability distribution or in other words, the Levy motion is transformed 
to the Brownian motion [||], ||. 

The possibility to develop the path integral over the paths of the Levy 
motion was discussed by Kac |7| , who pointed out that the Levy path integral 
generates the functional measure in the space of left (or right) continues 
functions having only discontinuities of the first kind. 

In the Refs.[§], it was shown that the path integral over the Levy- 
like quantum-mechanical paths allows to develop the generalization of the 
quantum mechanics. Namely, if the path integral over Brownian trajectories 
leads to the well known Schrodinger equation, then the path integral over 
Levy trajectories leads to the fractional Schrodinger equation. The fractional 
Schrodinger equation includes the space derivative of order a instead of the 
second (a = 2) order space derivative in the standard Schrodinger equation. 
Thus, the fractional Schrodinger equation is the fractional differential equa- 



tion in accordance with the modern terminology (see, for example, |T0|-||15||). 
This is the main point for the term fractional Schrodinger equation or for 
more general term fractional quantum mechanics, fQM ||. As it was men- 
tioned above at a = 2 the Levy motion becomes the Brownian motion. Thus, 
the fQM includes the standard QM as a particular Gaussian case at a = 2. 
Quantum mechanical path integral over the Levy paths at a = 2 becomes 
the well known Feynman path integral 0, [[J. 
The paper is organized as follows. 

In Sec. 2 the quantum mechanical path integral over the Levy paths has 
been discussed and the 3D fractional Schrodinger equation has been derived 
in term of the Riesz fractional derivative. 

In Sec. 3 we prove the hermiticity of fractional Hamilton operator in the 
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Hilbert space with scalar product defined by the same way as for conventional 
quantum mechanics. The parity conservation law for the fractional quantum 
mechanics has been established. 

Time-independent fractional Schrodinger equation has been introduced 
and its properties have been studied in Sec. 4. As physical applications of 
the time-independent fractional Schrodinger equation we have found (i) the 
energy spectrum and equation for the orbits radius for a hydrogen-like atom 
- fractional "Bohr atom", (ii) the energy spectrum of ID fractional oscillator 
in the semiclassical approximation. 

In Sec. 5 a new equation for the fractional probability current density has 
been derived and discussed. 

In the Conclusion we discuss the relationships between the fractional and 
the well-known Schrodinger equation. 



2 Path integral 

2.1 Path integral over the Levy paths 

If a particle at an initial time t a starts from the point r a and goes to a final 
point Tb at time tb, we will say simply that the particle goes from a to b and its 
path r(t) will have the property that r(t a ) = r a and r(tb) = r&. In quantum 
mechanics, then, we will have an quantum-mechanical amplitude, often called 
a kernel, to get from the point a to the point b. This will be the sum over all 
of the trajectories that go between that end points and of a contribution from 
each. If we have the quantum particle moving in the potential V(r) then the 
fractional quantum-mechanical amplitude KL(rbtb\r a t a ) may be written as 



K L {r b t b \r a t a ) = (1) 

r(t b )=r b ( . t b 

J Dr(r) J Dp(r) exp U / dr[p(r) r (r) - H a (p(r), r(r))] 

r(t a )=r a V ta 

where h is the Planck's constant, r denotes the time derivative, H a (p(r), r(r)) 
is the fractional Hamiltonian given by 

H a (p,r)=D a \p\ a + V(r), 1< a < 2, (2) 
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with the replacement p — > p(r), r — > r(r) and {p(t), r(r)} is the particle 
trajectory in phase space. The quantity D a has physical dimension 

[D a ] = erg 1 "" • cm" • sec"". 

r(t b )=r b 

The phase space path integral j Dr(r) /Dp(r)... in Eq.(P is defined 

r(t a )=r a 

by 



r(*b)=rb 

/ 

r(t a )=r a 



J Dr(r) j Dp(r)... = (3) 



oo oo 



Pi(ri-r Q ) .D Q ^|pi| a l [ ,pAr(r 6 - r Ar _i) .D a <;\pN 



a 



exp < z % > x...xexp < i i 



hi \ h h 



. . . . 



here ? = (t b - t a )/N. 

The exponential in Eq.([l]) can be written as exp{|S' Q ,(p, r)} if we in- 
troduce classical mechanical action for the trajectory {p(£), r(£)} in phase 
space 

t b 

S a (p, r) = J dr(p(r) r (r) - i7 a (p(<r), r(r), r). (4) 

When a = 2, D a = l/2m, (m is the mass of a particle) Eq.flU) is trans- 
formed into the well known Hamiltonian with the kinetic energy p 2 /2m and 
Eq.(fJ) becomes the definition of the Feymnam path integral in the phase 
space representation, see for example [T(d| . 



2.2 Fractional Schrodinger equation 

The kernel i£z, (r&tf,|r a i a ) which is defined by Eq. ([[]), describes the evolution 
of the quantum-mechanical system 
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■0/( r &> = J dr a K L (r b t b \r a t a ) ■ ?pi(r a , fa), (5) 

where ipi(r a ,t a ) is the wave function of initial (at the t = fa) state and 
ipf(ri,,tb) is the wave function of final (at the t = fa) state. 

In order to obtain the differential equation for the wave function ?/>(r, t) we 
apply Eq.fl5|) in the special case that the time differs only by an infinitesimal 
interval e from t a 

V>(r, t + e)= J dr'K L (r, t + e\r' , t) ■ i[>(r', t). 

t+e , 

Using Feynman's approximation / £ZrV(r(r),r) — eV( I ^ L ,t) and the 
definition given by Eq.([l]) we have 

if>(r,t + e) = 

n -p _|_ y f 

-j-eV(-^,t)}^(r>,t). 
hand sides in power series 



f , / 1 /• , r . P(r 7 -r) . A»e|p| a 

— oo 

We may expand the left-hand and the right- 



oo 



— oo 

(6) 

(l-^n^*))-v>(r',*). 

Then, taking into account the definitions of the Fourier transforms 
i>(r,t) = J dpe^(p(p,t), (p(p,t) = J dpe'^ip^t), 

and introducing the 3D quantum Riesz fractional derivative^ (—h 2 A) a ^ 2 



lr The Riesz fractional derivative was originally introduced in JT7| 
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(-h 2 A) a ^(r,t) = — J d 3 pe^\p\MP,t), (7) 
(here A = d 2 /dr 2 is the Laplacian) we obtain from Eq.(f|) 

^(r, t) + = V(r, t) - ^(-^ 2 A)^(r, t) - I e V(r, ftyfo t). 

This will be true to order e if ip(r,t) satisfies the differential equation 

lh ^A = D a (-h 2 A) a/2 ip(r, t) + V(r, t)^(r, t). (8) 

This is the fractional Schrodinger equation. The space derivative in this 
equation is of fractional (noninteger) order a. 

The above consideration is in fact the generalization of the well known 
Feynman approach to reduce the path integral to the differential equation 

1, !■ 

Equation (19) may be rewritten in the operator form, namely 

th^ = H a ij, (9) 
where H a is the fractional Hamiltonian operator 

H a = D a {-% 2 A) a ' 2 + V(v, t). (10) 

By definition ( |I0"D and introducing the momentum operator p = z^V one 
may obtain the fractional Hamiltonian H a in the form given by Eq.(@). 

Since the kernel KL(r b t b \r a t a ) thought of as a function of variables r b ,t b , 
is a special wave function (namely, that for a particle which starts at r a ,t a ), 
we see that Kl must also satisfy a fractional Schrodinger equation. Thus 
for the quantum system described by the fractional Hamiltonian Eq.([T0|) we 
have 



ih^-K L (r b t b \T a t a ) = D a (-h 2 A b ) a / 2 K L (r b t h \r a t a ) + V(r b , t)K L (r b t b \r a t a ), 
ot b 

(11) 

where t b > t a , and the low index "6" at the A b means that the fractional 
derivative acts on the variable r b . 
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3 Quantum Riesz fractional derivative 

3.1 Hermiticity of the fractional Hamilton operator 

The fractional Hamiltonian H a given by Eq.(|l(J) is the Hermitian operator 
in the space with scalar product 

(0, X )= / dT<f>*(r,t)x(r,t), (12) 

— oo 

where the sign * means as usual complex conjugate. 

To prove hermiticity of the fractional Hamilton H a let us note that in 
accordance with definition of the quantum Riesz fractional derivative given 
by Eq.flj]) there exists the integration by parts formula 

(0, (-h 2 A)^ x ) = ((-^A)^, x). (13) 
The average energy of fractional quantum system with Hamiltonian H a 



is 



E a = J drr(r,t)H a ^(r,t). (14) 

— oo 

Taking into account Eq. (|l3~D we have 

oo oo 

E a = J dn/>*(T,t)H a i/>(T,t) = J dr(H^(r,t)y^(r,t) = E* a . 

-co — oo 

and as a physical consequence, the energy of a system is real. Thus, the 
fractional Hamiltonian H a defined by Eq.(|T0|) is the Hermitian or self-adjoint 
operator in the space with the scalar product defined by Eq. ([12]) 



(H+<f>, X ) = & H aX ). (15) 

3.2 Parity conservation law for the fractional quantum 
mechanics 

It follows from the definition (0) of the quantum Riesz fractional derivative 
that 
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(-h 2 A) a ' 2 exp{z£p} = \p\ a exp{z^=}. (16) 

Thus, the function exp{ipx//i} is the eigenfunction of the 3D quantum 
Riesz fractional operator (— Ti 2 A) a l 2 with eigenvalue |p|°. 

The operator (— h 2 A) a ^ 2 is symmetrized fractional derivative, that is 

{-h 2 A r ) a/2 ... = {-h 2 A_ r ) a/2 .... (17) 

Because of the property ( O ) the fractional Hamiltonian H a given by 
Eq.(|T0|) remains invariant under inversion transformation. Inversion, or to 
be precise, spatial inversion consists in the simultaneous change in sign of all 
three spatial coordinates 

r — > — r, x — > —x, y — > —y, z — > —z. (18) 

Let us denote the inversion operator by P. The inverse symmetry is the 
fact that P and the fractional Hamiltonian H a commute, 

PH a = H a P. (19) 

We can divide the wave functions of quantum mechanical states with a 
well-defined eigenvalue of the operator P into two classes; (i) functions which 
are not changed when acted upon by the inversion operator, 

P^ + (r) = V+(r); 

the corresponding states are called even states; (ii) functions which change 
sign under the action of the inversion operator, 

Pip-(r) = —if)_(r)) 
the corresponding states are called odd states. 

Equation (|19|) express the "parity conservation law"; if the state of a 
closed fractional quantum mechanical system has a given parity (i.e. if it is 
even, or odd), then this parity is conserved. 
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4 Time-independent fractional Schrodinger equa- 
tion 



The special case when the Hamiltonian H a does not depend explicitly on 
the time is of great importance for physical applications. It is easily to see 
that in this case there exist the special solution of the fractional Schrodinger 
equation of the form 

^(r,t) = e-^0(r), (20) 

where 0(r) satisfies 

H a <f>(r) = E<f>(r), (21) 

or 

D a (-h 2 A) a/2 (j)(r) + V(r)(j)(r) = E(f>(r), 1< a < 2. (22) 

The equation (|22| ) we call by the time-independent (or stationary) frac- 
tional Schrodinger equation. 

4.1 Fractional Bohr atom 

The hydrogenlike potential energy V(r) is 




Then the fractional Schrodinger equation ( ]2"2] ) has a form 

Zp 2 

J D a (-^ 2 A)"/ 2 0(r) - — 0(r) = E<j>{v\ (23) 

and can be treated as fractional eigenvalue problem. 

The total energy of considered quantum mechanical system is 

E = E kin + V, 
where E kin is the kinetic energy 
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E kin = Da\p\ a , (24) 

and V is the potential energy 

y = m 

It is well known that if the potential energy is a homogeneous function 
of the co-ordinates and the motion takes place in a finite region of space, 
there exists a simple relation between the time average values of the kinetic 
and potential energies, known as the virial theorem (see, page 23, fl8f). It 
follows from the virial theorem that between average kinetic energy (|2"4] ) and 
average potential energy (^) there exist the relation 

aE km = -V, (26) 
where the average value / of any function of time is defined as 



oo 




In order to evaluate the energy spectrum of the fractional hydrogenlike 
atom let us remind the Niels Bohr postulates [I9j : 



1. The electron moves in orbits restricted by the requirement that the 
angular momentum be an integral multiple of h, that is, for circular orbits 
of radius a n , the electron momentum is restricted by 

pa n = nh, (7i = l,2,3,...), (27) 

and furthermore the electrons in these orbits do not radiate in spite of 
their acceleration. They were said to be in stationary states. 

2. Electrons can make discontinuous transitions from one allowed orbit 
corresponding to n = i%2 to another corresponding to n = ni, and the change 
in energy will appear as radiation with frequency 

oo = En2 ~ Eni , (na > n x ). (28) 

An atom may absorb radiation by having its electrons make a transition 
to a higher energy orbit. 
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Using the first Bohr's postulate and Eq.(p6|) yields 

fnh\ a Ze 2 

aD a — = , 

\a n J a n 

from which it follows the equation for the radius of the fractional Bohr 
orbits 

a n = a n^r, (29) 

here ao is the fractional Bohr radius (the radius of the lowest, n = 1 Bohr 
orbit) defined as, 

ao = \-^-r ) • (30) 
By using Eq.(p6[) we find for the total average energy E 

E = (1 - a)E km . 

Thus, for the energy levels of the fractional hydrogen-like atom we have 

E n = -(a - l)E n~^, 1 < a < 2, (31) 

where E is the binding energy of the electron in the lowest Bohr orbit, 
that is, the energy required to put it in a state with E = corresponding to 
n = oo, 



i 

a-1 



The energy (a — 1)E can be considered as generalization of the Rydberg 
constant of the standard quantum mechanics. It is easy to see that at a = 2 
the energy (a — 1)Eq is transformed into the well known expression for the 
Rydberg constant, Ry = me 4 /2h 2 . 

The frequency of the radiation uo associated with the transition, say, for 
example from k to n, k — > n, is, 



(a - l)E 



n < 



(k > n) (33) 



11 



The new equations (|29|)- (|33|) give generalization of the "Bohr atom" the- 
ory. In a special Gaussian case, a = 2 (standard quantum mechanics) 
Eqs.(P5[)-(P^) reproduce the well known results of the Bohr theory |TjJ, PD| ]. 
The existence of Eqs.(p^)-(|33|) is a result of deviation of fractal dimension 
^fractal °f the Levy-like quantum mechanical path from 2, d^^i = a < 2. 



4.2 Spectrum of the ID fractional oscillator in semi- 
classical approximation 

Fractional oscillator introduced in || is the model with the fractional Hamil- 
tonian operator H a /3 , 

H a ,p = D a (-h 2 A) a/2 + q 2 \rf, Ka<2, 1< < 2, (34) 

where r is the 3D vector, A = d 2 /dr 2 is the Laplacian, the operator 
(-h 2 A) a / 2 is defined by the Eq.(0) and q is a constant with physical dimen- 
sion [q] = erg 1 / 2 ■ cm"^ 2 . 

The ID fractional oscillator with the Hamilton function H a p = D a \p\ a + 
g 2 |x|^ poses an interesting problem for semiclassical treatment. We set the 
total energy equal to E, so that 

£ = A>r + g 2 |af , (35) 

whence 

l/a 



\P\ = {-^(E - q 2 \xf) 

At the turning points p = 0. Thus, classical motion is possible in the 
range \x\ < (E/q 2 ) 1 ^. 

A routine use of the Bohr-Sommerfeld quantization rule |22| yields 



2-Kh(n + -) = f pdx = 4 pdx = (E - q 2 \xf) 1/a dx, (36) 

o Da " o 

where the notation § means the integral over one complete period of the 
classical motion, x m = (E/q 2 ) 1 ^ 13 is the turning point of classical motion. To 
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evaluate the integral in the right hand of Eq.(|36|) we introduce a new variable 
y = x(E I ' q 2 ~)~ 1 ^ . Then we have 

J{E- q 2 \x\ P f /a dx = -^E^ J dy(l - y p ) 1/a . (37) 
o ^ o 

The integral over dy can be expressed in the terms of the S-functionQ 
Indeed, substitution z = y 13 yields 

J dy(l - /)V« = I j dzzl-\\ - z)$ = -B(- i + 1). (38) 



With help of Eqs.(|7|) and © we rewrite Eq.® as 

2nh(n + -) = —r^ E^\-B(\, - + 1). 

The above equation gives the value of the energies of stationary states for 
ID fractional oscillator, 

This new equation generalize the well known energy spectrum of the stan- 
dard quantum mechanical oscillator (see for example, ||22|| ) and is transformed 
to it at the special case a — 2, /3 — 2. 

As it follows from Eq . (p9|) at 

a(3 r (40) 



a + (3 

the energy spectrum becomes equidistant. When 1 < a < 2 and 1 < (3 < 2 
the condition given by Eq.(f40"D takes place for a = 2 and (3 = 2 only. It 
means that only standard quantum mechanical oscillator has the equidistant 
energy spectrum. 



The B{a 1 b) function has the familiar integral representation |2l| 
l 

B(a, b)= J duu^il - u) h ~ x . 
o 
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5 Current density 

By multiplying Eq.(||) from left by ip*(r,t) and the conjugate complex of 
Eq.(||) by ^(r, t) and subtracting the two resultant equations we finally obtain 

d 



Da 

ih 



-Jd 3 r(r(r } mr J t))= (41) 
J d 3 r (V»*(r, t)(-h 2 A) a/2 iP(r, t) - ^(r, t)(-h 2 A) a/2 ^*(r, t)) . 



From this integral relationship we are led to the following well known 
differential equation 

^M + di V j(r,t) = 0, (42) 

where 

p(r,t)=^MMr,t), (43) 

is the density of probability and the vector j(r, t) can be called by the 
fractional probability current density vector 



j(r, t) = ^ (V*(r, tX-tfA)"' 2 - 1 V^(r, t) - ^(r, t)(-% 2 A) a/2 ~ l Vip*(r, t) 

(44) 

where we use the following notation 

Introducing the momentum operator p = 4 V we can write the vector j 
in the form 

j =D a {ip(p 2 r /2 ~ l pr + v>*(p*t /2_1 p*v>) , i< « < 2. (45) 

When a = 2, D a = l/2m Eqs.(f0]), fl45|) becomes the well known equa- 
tions of the standard quantum mechanics (see, for example [f^Q. Thus we 
conclude that the new Eqs.([H]) and fl45|) are the fractional generalization of 
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the well known equations for probability current density vector of standard 
quantum mechanics. 

To this end, we express Eg. (f4"5D in the terms of the velocity operator, 
which is defined as usual 

where f is the operator of coordinate. Using the general quantum me- 
chanical rule for differentiation of operator 

d ~ i \xj i 
di T= h [H °> T] > 



we have 



v = j{H q y - TH a ). 
Further, with help of the equation 

df 

f(p)r-rf(p) = -iti-j^, 

which holds for any function f(p) of the momentum operator, and taking 
into account Eq. (0) for the Hamiltonian we find the equation for the velocity 
operator 

v = aD a \p 2 \ a/2 ~ l P, (46) 

here p is the momentum operator. By comparing of Eqs.(^) and ([ffi|) 
we finally conclude that 

j =- Uvil)* + tp*^) , 1 < a < 2. (47) 
a 

To get the probability current density equal 1 (the current when one 
particle pass through the unit area per unit time) the wave function of a free 
particle have to be normalized as follows 



^(r,t) = ^— exp{-pr--£t}, E = D a |p| a , 1< a < 2, (48) 
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where v is the particle velocity, v = aD a p a 1 . Indeed, by substituting 
Eq.(||) into Eq.© we find 

j=^, v = a J D Q |p 2 |f- 1 p, (49) 
that is the vector j is the unit vector. 

Equations (fHj), ( |45l) and (^7|) are the fractional generalization of the 
well known equations for probability current density vector of the standard 
quantum mechanics 



6 Conclusions 

Fractional generalization of the Schrodinger equation has been studied. We 
have established the integration by parts formula for the quantum Riesz frac- 
tional derivative and used it to prove hermiticity of the fractional Hamilton 
operator. The parity conservation law for the fractional quantum mechanics 
has been observed. The time-independent fractional Schrodinger equation 
has been introduced. As physical applications of the time-independent frac- 
tional Schrodinger equation we have found the energy spectrum and equation 
for the orbits radius for a hydrogen-like atom - fractional "Bohr atom". The 
energy spectrum of ID fractional oscillator has been obtained in the semi- 
classical approximation. 

The generalization of the fractional probability current density has been 
derived and discussed. 

The generalized equations (§, ©-(H), §§, §5|) and @7|) are 



transformed into the well known equations of conventional quantum mechan- 
ics if we put the Levy index a = 2. In other words the fractional quantum 
mechanics includes the standard quantum mechanics as a particular Gaus- 
sian case at a = 2. Quantum mechanical path integral over the Levy paths 
at a = 2 becomes the well known Feynman path integral and the fractional 
Schrodinger equation becomes the Schrodinger equation. 

The fractional Schrodinger equation provides us with general point of view 
on the relationship between statistical properties of quantum-mechanical 
path and structure of the fundamental equations of quantum mechanics. 
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